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Abstract 

We consider 5-dimensional spacetimes of constant 3-dimensional spatial cur- 
vature in the presence of a bulk cosmological constant. We find the general 
solution of such a configuration in the presence of a Gauss-Bonnet term. Two 
classes of non-trivial bulk solutions are found. The first class is valid only 
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under a fine tuning relation between the Gauss-Bonnet coupling constant and 
the cosmological constant of the bulk spacetime. The second class of solutions 
are static and are the extensions of the AdS-Schwarzchild black holes. Hence 
in the absence of a cosmological constant or if the fine tuning relation is not 



^ ' true, the generalised Birkhoff 's staticity theorem holds even in the presence of 



Gauss-Bonnet curvature terms. We examine the consequences in brane world 

cosmology obtaining the generalised Friedmann equations for a perfect fluid 

3-brane and discuss how this modifies the usual scenario. 
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I. INTRODUCTION 



The intriguing possibility that our Universe is only part of a higher dimensional spacetime 
JU, p| has raised a lot of interest in the physics community recently |J, ||, 0, [|7|], 
||. In particular 5 dimensional brane Universe models and their cosmology have been 
extensively studied (see for example 0, |K|, fllT] ). The Universe in this case is a gravitating 



homogeneous and isotropic brane or domain wall evolving in a constant negative curvature 
spacetime. 

Amongst the interesting features of such a toy-model configuration is that it verifies a 
generalised version of Birkhoff 's staticity theorem ||10|| , a constant curvature spacetime 



of constant 3-space curvature is locally static; more specifically an ADS-black hole solution 
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ds 2 = -h(r)dt 2 + h~\r)dr 2 + r 2 ( - X _ 2 + X 2 dn 2 H (1) 



dx 2 

1 - KX 2 

where h(r) = k — + k 2 r 2 , with k — 0, ±1 and /z, A; 2 are related to the black hole mass 
and bulk cosmological constant respectively. To ensure the validity of this theorem it is 
essential firstly that the brane Universe is of co-dimension 1 i.e. a domain wall type defect 
and secondly that the brane Universe is homogeneous and isotropic. The theorem in turn 
implies a certain number of physical properties for the configuration, in particular that 
the only dynamical degree of freedom is the wall's trajectory or equivalently, for the 4- 
dimensional observer stuck on the brane, the expansion rate of the Universe. Hence although 
we have introduced an extra dimension, the number of dynamical degrees of freedom does 
not alter with respect to standard 4 dimensional FLRW cosmology. Just like in 4 dimensional 
cosmology, once given an equation of state relating energy density and pressure one obtains 
the expansion rate or equivalently the brane Universe trajectory. It follows rather elegantly, 



TTJ, that it is totally equivalent to study a brane Universe evolving in a static background 
in the manner of [|K]] to a fixed brane Universe in a time-dependant background ||. 

When treating higher dimensional gravity theories we should keep in mind that 4 dimen- 
sional gravity is quite special for a numerous number of reasons. For example D = 4 gives 



the minimal number of dimensions where the graviton is non trivial and has exactly two 
polarisation degrees of freedom whereas at the same time gauge interactions of the Stan- 
dard Model are renormalisable. Another special property of 4 dimensional gravity is the 
uniqueness of the Einstein-Hilbert action. 

In D > 4, however, the situation is quite different. In 5 dimensions, in order to obtain the 
most general unique action, i.e. giving rise to a second order symmetric and divergence-free 
tensor, and to field equations that are of second order in the metric components, we have to 
add the Gauss-Bonnet term to the usual Einstein-Hilbert plus cosmological constant action. 
This is part of Lovelock's theorem ||14j| . Furthermore in an effective action approach of 
string theory, the Gauss-Bonnet term corresponds to the leading order quantum correction 



to gravity in particular in the case of the heterotic string fll5| . The Gauss-Bonnet coupling 
constant is related to the Regge slope parameter or string scale. Furthermore one of the 
important properties of string theories is that they contain no ghosts. Interestingly as was 
demonstrated in Jl6[ the only curvature squared terms to give ghost-free self-interactions for 
the graviton (around flat spacetime) is precisely the Gauss-Bonnet combination. 

The reason for all these nice properties shared by the Einstein-Hilbert and the Gauss- 
Bonnet terms can be understood from a purely geometrical point of view. The Gauss- 
Bonnet term is the generalised Euler characteristic of a 4 dimensional spacetime. It yields 
in D = 4 a boundary term hence a topological and not dynamical contribution. This is 
a quite general and elegant fact. Indeed we remind the reader that in a similar fashion 
the Einstein-Hilbert action in 2 dimensions is related to the usual Euler characteristic of 
a 2-dimensional manifold. For example in string field theory the Euler characteristic x is 
related to the string coupling constant g s , governing "surface diagrams" in the perturbative 
regime. In general every spacetime of even dimension 2n is accompanied by its generalised 
Euler characteristic; which we have to add to the gravitational action of a 2n + 1 manifold 
in order to preserve uniqueness. Thus for example in 10 dimensions one has the Euler 
characteristics of (cosmological constant), 2 (Einstein-Hilbert), 4 (Gauss- Bonnet), 6, and 
8 dimensional manifolds |14|]. So from this discussion it would seem natural to include the 



Gauss-Bonnet term in a 5 dimensional spacetime, all the more since we are interested in toy 
models merging string theory with standard cosmology 

Madore and collaborators have considered this term in order to stabilise the 5th dimen- 



sion in Kaluza-Klein theories JT7J whereas there was a lot of effort in the 80's to obtain exact 
solutions in Gauss-Bonnet theories in view to their relevance to quantum gravity corrections 
of string theory (see for instance [ 18] , |0| | 20fl , [p2| ). More recently in the context of 



brane universe it has been shown that the localised graviton zero mode persists in the RS 
model in the presence of a Gauss-Bonnet term f23j, |24j], [ffijfl . Cosmological consequences 



have also been studied in pq| . However, only particular solutions in the bulk have been 
considered. Here we shall attack the problem in its full generality. We shall first of all find 
and discuss the full bulk solutions, and then we shall investigate the brane cosmology they 
induce. Not surprisingly Birkhoff's theorem will be in the centre of our analysis and its 
physical consequences. 

In the next section we set up the basic ingredients of the problem. In Section III we 
solve by brute force the field equations and find the general solution for the bulk spacetime. 
In Section IV we discuss the relevance of the bulk solutions to brane Universe cosmology in 
5 dimensions. We conclude in section V. 

II. GENERAL SETTING 

Consider the following 5-dimensional action, 

M 3 r r n 

S=— J d 5 x [R + 12k 2 + a{R^ 5 R^ s - AR^R^ + R 2 )\ , (2) 

where M is the fundamental mass scale of the 5-dimensional theory, A = —6k 2 is the 
negative bulk cosmological constant and the Gauss-Bonnet coupling constant a of dimension 
(length) 2 ,which we leave free, is the additional physical parameter of the problem. Setting 
a = we obviously get the usual Einstein-Hilbert action with cosmological constant in 5 
dimensions. As we discussed in the introduction, just like the Einstein-Hilbert action with 
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cosmological constant is unique in 4 dimensions, the gravitational action (0) is unique in 5 
dimensions. To put it in a nutshell, (|2|) is the most general action that will yield second 
order partial differential equations with respect to the metric components in 5 dimensions. 
For this reason and for clarity we shall restrict ourselves to 5 dimensions. 

Let us now consider a spacetime with constant three-dimensional spatial curvature. A 
general metric can then be written, 

ds 2 = e 2l/ ^B(t, z)- 2 /\-dt 2 + dz 2 ) + B(t, zfl z i-^— 2 + X 2 dnjj) (3) 

where B(t, z) and u(t, z) are the unknown component fields of the metric and k — 0, ±1 is the 
normalised curvature of the 3-dimensional homogeneous and isotropic surfaces. We choose 
to use the conformal gauge in order to take full advantage of the 2-dimensional conformal 
transformations in the t — z plane. This is the setup for a cosmological wall or brane-Universe 
of co-dimension 1. We note on passing that a co-dimension 2 or higher set-up would have 
lost the two dimensional conformal freedom. We will see in what follows that this freedom 
is essential for the integrability of the system wether or not we include the Gauss-Bonnet 
combination. 

The field equations we are seeking to solve are found by varying the above action (^) 
with respect to the background metric and read 



E 



Guv - §k 2 g^ v - a 



^(R a ^ s R^ s - AR lS W s + R 2 ) 



— 2RR^ V + 4R^R\ + AR^sR 1 ^ — 2R^s\RJ f6x = (4) 

where now the symmetric tensor e^ has replaced the usual Einstein tensor and is also 
divergence free. Taking the trace of (f|) one can show that for a solution, the Ricci scalar is 
a multiple of the Lagrangian in ([J) (see for instance ||21|| ). Thus the behaviour (in particular 
singularities) of the scalar curvature R is shared by the Gauss-Bonnet scalar in (|2]). Hence 
not surprisingly we can deduce that although the field equations change radically, spacetime 
curvature still plays the same physical role for the critical points of the action (Q). 
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III. THE GENERAL SOLUTION FOR THE BULK SPACETIME 



Before plunging in the field equations^ it is rather useful to review the generalisation of 



Birkhoff's theorem in the presence of a cosmological constant as it appeared recently in [3 I 
Furthermore we shall use exactly the same method to derive the general solution. 
In this subcase the field equations are obtained setting a = in (0) and read 

_2A 

There are two key ingredients in this method. First of all in order to make use of the t — z 
conformal symmetries of @ it is important to pass to light-cone coordinates, 

t — z t + z 

Secondly taking the combination R tt +R zz ±2R tz = 0, one obtains the integrability conditions 
which read[| 

B, uu - 2v u B jU = 0, (6) 
B >vv - 2u v B ;V = 0. (7) 

Note then that these are ordinary differential equations with respect to u and v respectively 
and are independent of the physical parameter of the problem, namely, the cosmological 
constant A. As their name indicates they are directly integrable giving 

B = B(U + V) e 2u = B'U'V (8) 

where U = U(u) and V = V(v) are arbitrary functions of u and v, and a prime stands for 
the total derivative of the function with respect to its unique variable. Using a conformal 
t r ansf or mat ion , 



1 For the full field equations see Appendix. 

2 From now on B :U represents the partial derivative of B with respect to u etc. 
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gives that the solution is locally static B = B(z) and Birkhoff's theorem is therefore true. 
Starting from a general time and space dependant metric, spacetime has been shown to be 
locally static or equivelantly that there exists a timelike Killing vector field (here Jj). Note 
that we did not have to find the precise form of the solution for B. The integrability condi- 
tions actually suffice to prove staticity and thus Birkhoff's theorem. By use of the remaining 
field equations we can then find the form of B, leading after coordinate transformation to 
the topological black hole solution ([!]). Note that the solution becomes t-dependent as we 
cross the event horizon of the black hole. For more details the reader can consult ||T] • 



Let us now turn to our case of interest with a ^0. In the presence of the Gauss-Bonnet 
term we can expect that if the system is indeed integrable then some integrability equation 
should be reproduced. Putting away technicalities this is the essence of what we shall do 
here. In analogy to the previous case let us take the combination, e u + e zz ± 2e <2 = 0. On 
passing to light cone coordinates @ we get after some manipulations 

(9B 4/3 e 2u + 36aKB 2/3 e 2v + AaB jU B jV ) (B, au - 2u u B >u ) = 

(9B 4/3 e 2u + 36aKB 2/3 e 2u + AaB >u B tV ) (B >vv - 2u v B >v ) = (9) 

Note how the Gauss-Bonnet terms factorise nicely leaving the integrability equations (|5|) we 
had in the absence of a. 

Let us neglect for the moment the degenerate case where either B u = or B v = 
corresponding to flat solutions [^7] (see Appendix). For B jU ^ and B v ^ the situation 
is clear: either we have static solutions and Birkhoff's theorem holds as in the case above or 
we will have 

*> = Aa ^~ B V (10) 

9^/3(52/3 + 4^) ^ V > 

Let us first examine the latter case, that we will call Class I solution. The two remaining 
field equations e xx = and e u — e zz = give after some algebra the simple relation, 

8ak 2 = 1 (11) 
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This is quite remarkable: if the coupling constants obey this simple relation (|TT|) then the 



B field is an arbitrary function of space and time. Note in passing that Class I solutions 
exist in arbitrary dimension d if the fine tuning relation (^jfw|~2) = 1 is satisfied. We can 
already deduce that Birkhoff's theorem does not hold for non zero cosmological constant. 
However in the absence of a cosmological constant it is always trivially true since (JTTJ) is 



impossible. Also we can note from ([LI]) that a positive Gauss-Bonnet constant a > 0, as 
in heterotic string theory, demands a negative cosmological constant and vice-versa. The 
Class I metric reads, 

= 9 ^fL) ^ + <*> + 52/3 (A + (i2) 

under the constraint (p]) where we emphasize that B(t, z) is an arbitrary function of t and 
z. To simplify somewhat set B = R 3 to get, 



<fa2 = R l jS (_^2 + dz 2 } + ^2 + ) (13) 



This solution has generically a curvature singularity for R z = ±-R 4 . The parameter a is 
related here to the 5-dimensional cosmological constant via ([TT]). The Class I static solutions 
are given by, 

ds ^-^ ie+ ^ +R2 {A + ^) (i4 » 

K +4^ K +4^ V 1 ^^ / 

with A = A(R) now an arbitrary function of R. Time-dependent solutions for a > are 
only possible for R 2 < 4a and k = — 1.0 

In order to obtain t and z dependent solutions it suffices to take the functional R to be 
a non-harmonic function. Take for instance R = exp(f(t) + g{z)), with / and g arbitrary 
functions of a timelike and spacelike coordinate respectively. Let us also assume k = for 
simplicity, the Class I metric in proper time reads, 



3 Note however that for a non-zero charge Q and spherical symmetry (k = 1) Birkhoff's theorem 

is always true as was shown by Wiltshire (see also pCj] ) 
4 For a < the situation is interchanged with static solutions possible only for k = 1 and R 2 < 

-4a. 



^-^t^- 2(/+9, (t^H- (15) 

Note here again that / is an arbitrary function of time. 

On the other hand if (|TI|) does not hold then Birkhoff's theorem remains true in the 
presence of the Gauss-Bonnet terms i.e. the general solution assuming the presence of a 
cosmological constant in the bulk and 3 dimensional constant curvature surfaces is static if 
and only if ( f71| j is not satisfied. In this case the remaining two equations give the same 
ordinary differential equation for B(U + V) which after one integration reads, 

B' + 9B 2 ^(k 2 B 2 ^ + K) + 18a(^^ + ^j = 9/i (16) 

where fi is an arbitrary integration constant. Then by making B the spatial coordinate and 
setting _B 1//3 = r we get the solution discovered and discussed in detail by Boulware-Deser 
k = 1) and Cai || (k = 0, -1),Q 



df = -V(r)df + *1 + r< i & + X'dah) (17) 



where V{r) — k, + |^[1 ± yjl — Sak 2 + 8^f], and /i plays the role of the gravitational mass. 
The maximally symmetric solutions are obtained by setting \i = 0. There are two AdS 
branches permitted by the solution for a > ( ||28|| , ||L8]1 ). We can have both de-Sitter and 
Anti-de-Sitter for a < 0. Generically, as shown by Boulware and Deser ||18|| , only one of the 
branches is physical, the '+' branch being classically unstable to small perturbations and 
yielding a graviton ghost. For a > and the '— ' branch there is a black hole singularity at 
r = 0, a unique event horizon and asymptotically one approaches the 5-d topological black 



hole solutions [12], pi (see also |p0||). For a general and thorough analysis of the Gauss 



Bonnet black hole solutions, and their thermodynamics we refer the reader to [ 28| , [18| 0. It is 



interesting to point out that only the planar black hole k = obeys the entropy-area formula. 



5 We have kept the same label as in (||) for the rescaled time coordinate . 

6 Perturbative AdS black holes with R 2 curvature terms and their thermodynamics have been 



discussed in |3lj| ) 
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Indeed it turns out, that the planar black hole shares exactly the same thermodynamic^] 
properties as the planar topological black hole (a = 0) although the two solutions differ 
considerably. We will come back to this point in the next section. Furthermore for small a 
we have, V(r) = « + k 2 r 2 (l + 2ak 2 ) - #[1 + 2a(2k 2 - £)] + 0(a 2 ) and indeed for a = we 
get the usual Kottler solution |p"2| . 

Now notice how (|TTD is a particular 'end' point for fll7|) since the maximally symmetric 
solution is defined only for 1 > 8ak 2 (for a < there is no such restriction). We can deduce 
in all generality that for 1 > 8ak 2 there is a unique static solution (|17D . When ( |TT1) is 

2 

satisfied and /i = 0, the two branches coincide and V = k + |-, which is then a particular 
Class I solution ([Tj|) for the value A(r) = V(r). For 1 < 8ak 2 no solutions exist. 



IV. BRANE WORLD COSMOLOGY 

Having evaluated the general solution in the bulk we now consider a 4-dimensional 3- 
brane where matter is confined. We furthermore suppose following the symmetries of our 
metric (f|) that matter on the brane is modelled by a perfect fluid of energy density p and 
pressure p. The brane is fixed at z = 0, and the energy-momentum tensor associated with 
the brane takes the form, 

T; (6) = ^idiag(-p(t),p(t),p(t),p(t),0). 

and the field equations read, 

v = m- 3 t( 6 >(4 (is) 

We will assume Z2 symmetry across the location of the brane at z — 0, and set M 3 = 1 for 
the time being. 

Now before proceeding there are three important points to take into account. First of all 
the Israel junction conditions are no longer valid since we have included the Gauss-Bonnet 

7 The interested reader could consult the more general approach of (2^] 
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term in the gravitational action. Although the Gauss-Codazzi integrability conditions are 
universal for any spacelike or timelike hypersurface (see for instance ||32||), the Israel junc- 
tion conditions have to be generalised in order to take into account the addition of the 
Gauss-Bonnet term |33| in the gravitational action (FJ). So in order to evaluate the brane 
junction conditions we choose to bifurcate the Israel junction conditions, integrating the 
field equations (|18|) on an infinitely small interval across the brane location at z = 0. 

The second important point is that since the field equations are of second order we will 
encounter at most second derivatives of z and therefore the metric component fields have to 
be continuous. Indeed first order derivatives contain a jump in the metric given by means 
of the Heaviside distribution whereas second order derivatives contain a Dirac distribution 
at z = 0, to be matched with the brane energy- momentum tensor (0). Note that had we 
considered any other combination of quadratic curvature terms in the action, the situation 
would have been different. The good behavior of the Gauss-Bonnet combination is coherent 
with the fact that (fj) is unique in 5 dimensions, just as ordinary Einstein Hilbert gravity 
plus cosmological constant is unique in 4 dimensions. Hence we can expect a regular gravity 
theory and hence regular boundary conditions. 

The final remark turns out to be crucial for the correct evaluation of the junction condi- 



tions and has been a source of confusion in the related literature. Indeed note that, in (|T8|). 
first order derivatives with respect to z, multiplying second order derivatives of the metric 
functions, are always encountered as squares. Thus, the first order part involving Heaviside 
distributions turns out to be equal to +1 everywhere except at z = 0. Although this is a 
removable discontinuity it occurs in the location of the Dirac distribution and hence the 
junction conditions are not obtained by simply matching the Dirac distributions in the field 
equations It is imperative that we integrate over an infinitely small interval across the 
brane location Q and then take the limit. In doing so the t — t and x~X components of ([18]) 
give respectively the energy density and pressure on the brane: 



8 We thank Stephen Davis for discussions on this point 
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P(t) = - 2 -e^B-*BMi + yBj) (19) 



p(t) = \e^B-- 



h(Bu z - ~B >Z ) + B !Z [6e 2u B^ 3 + SaB(B >tt - u t B >t ) - ^B 2 t ] 



(20) 



with Ji = gS^e 2 " + 36a/t.B 2/3 e 2iy - 4a(B 2 a - B 2 t ) (see appendix). All the functions in the 
RHS of (|19"D and (^0|) are evaluated at z = + since we have assumed ^-symmetry. The 
domain wall case (p = —p) corresponding to a Poincare invariant brane has been treated in 



24 , 25 



We will first focus on the Class II solution in the bulk. In this case the junction conditions 
remain invariant under the conformal boost u — > f(u), v — > f(v) just as for a = 0. Therefore 
for a fixed brane (or boundary) at z = there is a single degree of freedom U' or V for the 
bulk spacetime which is evolving in time. By virtue of Birkhoff 's theorem this is equivelant 
to taking a moving brane (or boundary) in the static black hole background. We make use 
of this fact now to pass on to the static bulk configuration (for a detailed discussion see 

(0))- 

Consider a brane Universe observer. The expansion parameter (or wall's trajectory) 
reads, R(t) = S 1 / 3 (t, 0) whereas proper time r is given by dr = e u ^B- 1 / 3 (t,0)dt. For the 
solution (|I7p we have relations @ and for example the Hubble expansion rate is given by, 

_ IdR _(U' + V')B' 

First, using (|8|), (16), ( |19D and ( p0|) , we may obtain the standard conservation equation on 



the brane: 

^ + 3H(p + p) = (21) 

which is a consequence of the Bianchi and Bach-Lanczos identities for the Einstein tensor 
and Gauss-Bonnet terms respectively in (ffl). Then, using (H), flTED and (IB5p we get the 
generalised Friedmann equation: 

{&H*+Th<>{*+T)'M* + ™) (22) 
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where we have defined 



and V(R) is the black hole potential (0). This equation relates the brane trajectory R(t) 



with the energy density of the brane. Equations (22) and (^l|) fix the unique degree of 
freedom in the bulk R = R(t). In Einstein-Hilbert brane cosmology where a = 0, p 2 depends 
linearly on H 2 + ^r, where h(R) is the Kottler potential ([!]). Indeed taking heuristicaly the 
limit of small Gauss-Bonnet coupling (a — > 0) in ( p2|) yields the usual Friedmann equation 
for a 3-brane embeded into a (Einstein-Hilbert) five-dimensional bulk (see for instance |fL0|j): 

H 2 = ^- — -k 2 + -^- + 0(a) (23) 
36 R 2 R 4 K J K ' 

for the lower '+' sign in (^). This sign corresponds to the stable branch for the black hole 
solution (|17D as demonstrated by Boulware and Deser f\ At early times, ( p3f ) leads to 
unconventional cosmology H 2 oc p 2 0, in contrast with standard four-dimensional cosmol- 
ogy where H 2 oc p. Now, note from ( p2|) how Gauss-Bonnet gravity yields p 2 depending 
also on (H 2 + ^r^) 2 and (H 2 + ^P^) 3 - In general the three powers in the RHS of (|22|) are 
expected to dominate successively during the cosmological evolution of the universe. This 
may have interesting consequences for early as for late time brane cosmology. Generically, 
the cosmological evolution resulting from (|22|) strongly depends on the epoch under con- 
sideration and on the order of magnitude of the bulk lagragian parameters a and k. A 
detailed study lies beyond the scope of a simple application of Birkhoff 's theorem to higher 
dimensional theories and will be undertaken in future work. 

As an illustrative example we consider here the late time cosmology of a spatially flat 
(k = 0) expanding universe, for the particular relation between the bulk parameters 

Sak 2 = -3 , (24) 



9 Note that for a — > the upper '+' branch of ([l?]) yields a singular negative term in the RHS of 
|2^), which ties in nicely with the results of ]l8| showing the instability of this branch. 
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We will consider the '— ' (upper) sign in (^) and expand up to 0(-^) (large scale factor). 
In this particular case, both the effective cosmological constant on the brane and the black 
hole term Q in -j^ vanish. Equation ( p2|) then reduces to 

-aH 6 + H A - -^-H 2 = -^ + 0(4) 
3 2R A 192a K R 8 ' 

At sufficiently late times, the leading contribution for the brane energy density is if 2 oc p, as 
in standard four- dimensional cosmology, without the need to introduce any brane tension. 
For a > we have a positive cosmological constant in the five-dimensional bulk action, but 
however the bulk space-time is AdS, as may be seen through the expression of the black- 
hole potential V(R) fll7|) . Domain wall solutions (p = —p) without brane tension and in the 
presence of higher curvature terms, have been studied in |25| . It has been proved that they 



may allow for a massless normalisable four-dimensional graviton. We can regard fl24|) as a 
special relation involving the bulk parameters in contrast to the Randall- Sundrum relation 
relating brane and bulk parameters. The catch however is that we used the '+' (upper) 
branch of the solution in (|T7D, which is unstable according to ]Tj|. It would be interesting to 
study the cosmology of similar cases in higher dimensions with higher order Euler densities, 
which may be stable P5 . 



Let us now consider the case 8ak 2 = 1 and p = 0. Then, ( 22j) becomes: 



H 2 = (yp) 2/3 - 2k 2 - £ 

This is the result we obtain for a Class I solution in the bulk, by directly using (pTOj ) and 
flnj|) (extended for k ^ 0) in (|I9|). In this case, with (|20|), one also obtains (|21|). Again for 
Class I solutions, there is only one dynamical degree of freedom in brane cosmology, namely 
the function f(r) appearing in (p~5|) . 



10 This term is usually referred to as the 'dark radiation term' due its 'radiation' like behavior. 
This is however misleading since the bulk solution does not radiate (Birkhoff 's theorem) and there 
are furthermore no 'radiation' like particles in the bulk. 
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Now it is interesting to study the late time cosmology of an expanding brane-universe 
resulting from the generalised Friedmann equation (p2|). In doing so, we take into account 
the tension (vacuum energy) T of the brane (p —>■ T + p), and keep only linear terms in p 
and (large scale factor). We concentrate on the physical Class II solutions ('-' branch in 
(|T7|)), with a negative cosmo logical constant in the bulk action and with a > (as required 
by string theory), so that < 8ak 2 << 1. For a zero effective cosmological constant on the 
brane, one has to impose the modified Randall- Sundrum fine-tuning condition: 



1 - VI - 8ak 2 



[2 + Vl- Sak 2 ) (25) 



a 

which indeed allows for a Kaluza-Klein zero-mode localized on the brane and a finite volume 
element, rripi < oo |2~5| . Then, up to O(p) and O(j^), ( p2[ ) gives 



H 2 = ^f P - £ + (2 + VT^F)- 1 ^ (26) 
where the reduced four-dimensional Planck mass is given by 



and we have restored the fundamental mass scale M of the 5-dimensional theory. Note that 
the "cosmological" Planck mass, as defined above, agrees with the four-dimensional Newton 
constant obtained through estimation of the static gravitational potential at long distances 



along the brane [24 



Repeating the same procedure as above in the special case 8a k 2 = 1, one finds that ( p5|) 
and ( PT[ ) still hold, while the ^t-dependent term in the Friedmann equation fl2"6j ) turns out 



to be in Thus for < 8ak 2 < 1, one sees from ( |2"7| ) that: 



M 3 < m 2 Pl k 



whereas strict equality holds in the absence of the Gauss-Bonnet term (which is the leading 
quantum gravity correction term). Hence for fixed 4-dimensional Planck mass and cosmo- 
logical constant in the bulk, "quantum corrections" for gravity in the bulk tend to decrease 
the fundamental mass scale M of the 5-dimensional theory. 
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V. CONCLUSIONS 



In this paper we have studied Gauss-Bonnet brane cosmology in a 5 dimensional space- 
time. Our main motivation for including the Gauss-Bonnet term is that the usual 5 di- 
mensional gravitational action (|2]) is then unique [13] as we noted in the Introduction. 
Furthermore the Gauss-Bonnet coupling constant a provides a window to the leading quan- 
tum gravity correction coming from string theory. Throughout our analysis the technical 
difficulties induced by the inclusion of the higher order curvature term where seen to be 
overcome quite elegantly for the Gauss-Bonnet combination. 

Indeed starting from a homogeneous and isotropic 3-space in constant bulk curvature 
we found the general spacetime solutions to the field equations. Under a particular relation 
between the bulk cosmological constant and the Gauss-Bonnet coupling, a space and time 
dependant solution (Class I) of the field equations was found. If however this special relation 
is not satisfied then the unique solution is the black hole solution discovered and discussed 



m 



18], [21], [19] and [28|. Therefore quite elegantly Birkhoff's staticity theorem holds and 



all its interesting properties go through just like in the ordinary a = case [11 



As a concrete application to the generalised Birkhoff's theorem we studied brane cos- 
mology in 5 dimensions for a 4-dimensional perfect-fluid brane. As it turns out special care 
has to be taken when deriving the generalised Friedmann equations for the brane. Although 
there are no ill defined distributional products (unlike any other higher order curvature 
theory) a limiting procedure has to be undertaken in order to obtain the junction condi- 
tions. On doing so it is found that the generalised Friedmann equation involves a third 
order polynomial in H 2 which yields drastic changes to conventional (a = 0) brane cosmol- 
ogy. Also generically Gauss-Bonnet gravity tends to decrease the 5 dimensional fundamental 
mass scale, which is interesting if we interpret the Gauss Bonnet term as the leading string 
quantum gravity correction. It is now important to investigate whether the Gauss-Bonnet 
term can give ordinary late time FLRW cosmology without the usual fine tuning conditions 
needed in Einstein-Hilbert brane cosmology. Work in this direction is under way. 



16 



ACKNOWLEDGEMENTS 



It is a great pleasure to thank Brandon Carter, Stephen Davis, Emilian Dudas and John 
Madore for encouraging and discouraging remarks in the early stages of this work. We 
thank in particular Pierre Binetruy and Jihad Mourad for fruitful discussions. CC also 
thanks Ruth Gregory for pointing out [|2?J and Simon Ross. 

Note added: [Q has been brought to our attention, which appeared quasi-simultaneously 
to our paper and discusses some related issues. 

APPENDIX: FIELD EQUATIONS 



The field equations obtained from read, 

E xx = -g^Tg (B,u ~ B, zz )[l2e 2 »B 2 ^(2aK + B 2 ' 3 ) - h] 

+ ^-[B iZZ B itt + {yl - v 2 z ){B 2 t - B 2 Z ) + B >tz {2v tt B >z + 2B )t u z - B >tz ) 



9e 

- (B itt + B >zz )(p t B tt + v z B iZ )) 

20 « /r?2 D 2N2 h 



81e^B 2X '* > zJ 9e 4v 
~ 6 ^ 2/3 " " 6e2 ^ 2/3 ( 6aK + S2/3 )] = ( Al ) 

1 1 . 12e 2u k 2 

E tt ~ E zz - ^ t 2v B llZ ( B > U ~ B ' zz ) + £2/3 

6ne 2u _ 8a 2 _ 2 2 

54/3 27B w / 3 e 2u '* ' 

- l^rM - #,) = (A2) 
The integrability conditions e u + e zz ± 2s tz are: 

h {B, tt + B. zz + 2B M - 2u t B tt - 2v z B, z + 2v t B, z + 2B >t u z ) = (A3) 

h(B )U + B >zz - 2B >tz - 2v >t B it - 2p >z B >z - 2v±B. z - 2B, t v z ) = (A4) 

where h = §B^e 2v + 3QanB 2 ^ 3 e 2u - Aa{B 2 z - B 2 ). 
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In the degenerate case where either B jU = or B tV = (Class I solutions according to 



Taub) there are now two subcases. Either we obtain the Class I solution of Taub ]27[ which 
is simply flat Minkowski spacetime or we obtain, 



ds 2 = e^i-Aany^-dt 2 + dz 2 ) + (-4a«) + X 2 ^) (A5) 



1 - KX 2 

under once again (|TT|) . Note once more that v(t, z) is an arbitrary function of t and z and 
planar symmetry is not permitted. 
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